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Abstract 

We here investigate a cosmological model in which three fluids interact with each other 
involving certain coupling parameters and energy exchange rates. The motivation of the 
problem stems from the puzzling 'triple coincidence problem' which naively asks why the 
cosmic energy densities of matter, radiation and dark energy are almost of the same order 
of magnitude at the present time. In our model, we determine the conditions under triple 
interacting fluids will cross the phantom divide. 
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1 Introduction 



Despite several successes of the standard big bang cosmology based on Friedmann- Robert son- 
Walker (FRW) model, still a series of problems to be resolved like the horizon problem, 
flatness problem, dark matter (or missing mass) problem, structure formation, topologi- 
cal defects, matter-antimatter asymmetry and the cosmic coincidence problem etc. Most 
of these are related with the cosmic past of the observable universe while the cosmic 
coincidence problem has its origin in the recent time since it naively asks why certain cos- 
mological phenomena are occurring in our presence or in our times. Recent astrophysical 
observations give convincing evidence of an accelerating universe caused by dark energy 
characterized by the equation of state (EoS) parameter a; ~ — 1. It is yet unknown why 
the present energy density of the dark energy is approximately equal to that of dust-like 
matter. It is termed as the cosmic coincidence problem [1] . Till now several models have 
been proposed in an attempt to solve this problem such as the 'tracker field' [2], oscillat- 
ing dark energy |3] and the variable constants approach to name a few. It appears 
that the energy density of the radiation component is also almost equivalent to that of 
the matter and the dark energy i.e. Pm ~ Pr ~ Pa ot flm ^ fir ^ ^a, the so-called 
'cosmic triple coincidence problem' [5]. The question is: why this happens in the current 
or in recent times. The history of the parameter u suggests that it is no more a constant 
but possesses a parametric form uj{z), where z is the redshift parameter. Thus in the 
past uj = 1/3 corresponds to radiation and then uj = for matter. Later it evolved to 
quintessence < — 1/3 to cosmological constant uj = —1. This behavior suggests that in 
future, CO will be super- negative i.e. u < —1, which corresponds to the phantom energy. 
Thus in totality, we have the transition from a; > — 1 to < — 1 the so-called 'phantom 
crossing' or 'phantom divide' scenario, while we are observing = — 1 at the current time 
[6]. The coincidence problem in this context is rephrased as 'why is c<j = — 1 now?' 

In recent years, the usual coincidence problem is addressed by proposing an exotic 
interaction between dark energy and matter in which energy from pA is diluted or decayed 
into the p^ [3 El El [ini [ni [El [131 [11 [151 [ini [HI [1^ It is recently proposed that if these 
two components interact then some energy might dissipate into a third component px 
which is as yet hypothetical [T9]. The third component can be known form of matter 
or an altogether exotic fluid in which case some new physics will be required to explain 
the interaction. If we assume p^ = pr then the interaction between three fluids i.e. 
matter, radiation and dark energy will be quite interesting. It is well-known that matter 
and radiation were decoupled at the time of emission of cosmic microwave background 
(CMB) radiation at a redshift z ~ 1100. Thus both matter and radiation are almost non- 
interacting components but it can be anticipated that these two components do supposedly 
interact with the dark energy. Thus if dark energy and matter interact, the energy 
dissipated in the interaction is assumed to transfer to the radiation component and vice 
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versa for the radiation and the dark energy interaction. This dynamic interaction than 
hugely alters the effective equation of state of the three interacting fluids. Our analysis in 
this paper suggests that effective EoS of the three fluids will be interlinked to each other 
considerably and also constrained by coupling parameters. Moreover, the interaction 
naturally leads to the phantom crossing scenario. 

We here present a model in which three cosmic fluids (radiation, matter and dark energy) 
interact with each other with equal degree of freedom. Hence there are three coupling 
parameters involved which can take arbitrary positive or negative values not all zero 
simultaneously. The positive and negative values correspond to back and forth nature of 
energy exchange between the fluids. We emphasis here that the coupling parameters in 
our model are only constrained by the choices of effective EoS of the fluids. Thus in our 
model, exact EoS of the fluids hke in the non-interacting FRW model is not possible. As 
the exact EoS for the dark energy is unknown, its interaction with other fluids creates 
ambiguities in the determination of the exact EoS of the other fluids. Hence we also stress 
that precise values of the EoS's can be deduced only empirically from the phenomenology 
of the interacting fluids. 



2 The model of triple interacting fluids 

We start by assuming the background to be spatially flat, homogeneous and isotropic 
FRW spacetime 

ds^ = -df + a^it) [dr^ + r\de^ + sin^ ed(f)^)] , (1) 

where a{t) is the scale factor. We consider three fluids having the equations of state 
(EoS) Pi — uJiPii i — 1,2,3 where pi and pi are the corresponding pressures and the energy 
densities of the fluids, respectively. Also uji arc the dimcnsionlcss EoS parameters. For 
the sake of simplicity, we assume the three fluids to be perfect fluid-like since in general, 
interaction between fluids might lead to local inhomogeneities, which are ignored in the 
present paper. The equations governing the interaction of three fluids are expressed as 

Pi + 3H(l + u;i)pi = Q3-Q2, (2) 
P2 + 3H(l + u;2)p2 = Q1-Q3, (3) 
P3 + 3H{l + u;3)p3 = Q2-Q1. (4) 

Here Qi are the energy exchange (or dissipative) terms to be put ad hoc in the above 
equations. The explicit form of Qi should be determined only from the phenomenological 
and empirical results. However, from dimensional considerations, the quantity Qi should 
have dimensions of density into the time inverse. Choosing the later one to be Hubble 
parameter, we notice that Qi can be of the following forms: Qi ~ Hpi. This approximation 



3 



can be saturated to equality by inserting a dimensionless parameter (say Aj). Thus we 
can write ^20j (see [21] for more exotic expressions for Qi) 

= A,i/pioa-3(i+-«), 2 = 1,2,3. (5) 

In the last expression, Aj are the coupling constants which can take positive or negative 
values to yield two-sided energy exchange rather than one-sided. Also Eqs. (2) to (5) 
show that this a coupled system of three differential equations which needs to be solved. 
Further pio are the constant energy densities at some reference time t = to- Also H = a/a 
is Hubble parameter which determines the rate of expansion of the universe. Sum of Eqs. 
(2) to (4) yield the combined energy conservation 

3 

Y,[pi + 3H{l + uj,)p,] = 0, (6) 

i=l 

Also the energy conservation for the individual component (for the case of non-interacting 
fluids) yields 

p[ = p^oa-'^'+''^\ z = l,2,3. (7) 

Here pio are integration constants. Combining Eqs. (2), (3) and (4), we arrive at the 
density evolution of the interacting fluids as 

3 UJ2 — uji 3 uji — Us 

3 UJ-i — UJ2 3 UJ2 — UJi 

,3 . c.,-^a...^h^^:^^^^:^, (10) 

3 cji — 3 0^3 — u;2 

Here Cj's are constants of integration. Now addition of Eqs. (8), (9) and (10) results in 

p = p[X,+p',X2 + p',X3, (11) 

where 

X, = 1 + ^^ ^ (12) 

X. = 1 + ^. ^ (13) 

X, = 1 + ^^ ^ (14) 

3 [uJi - UJsj[uJ3 - UJ2) 

Here we have assumed Cj = pio and p = pi + p2 + Ps is the total energy density of the 
interacting fluids. Note that the case of non-interacting fluids is obtained by choosing 
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Ai = A2 = As = which yield Xi = X2 = = 1 i.e. p = p'l + P2 + Ps- Using the first 
FRW equation 

= -^p. (15) 

Differentiating Eq. (7) w.r.t t, we obtain 

p'. = -3{l + uJi)Hp,oa-'^'-^''^\ 2 = 1,2,3. (16) 

Differentiating Eq. (15) w.r.t t and then using Eqs. (11) and (16), we find 

-1/2 



H = -HV24^j2[p'Al+u;,)xJj2p:x) 

1=1 \i=l / 



(17) 



2.1 Role of parameter H 

The parameter H is significant as its possible signature governs the dynamics of the 
universe |3T]. The slowing down in the expansion takes place when 

^ 1 \ ^ 3(a;i -a;3)(a;2 -^1) 
uji > -1, Ai > , (18) 

LJ2 > -1, A2 > , (19) 

Ui - UJ3 

^ T X ^ 3(^1 -^3) (^3 -^2) 

UJ3 > -1, A3 > . (20) 

UJ2 - UJi 

Thus H will be negative when all the Ui > —1. This result corresponds to the quintessence 
dominated universe. 

Note that the vanishing H in Eq. (17) will yield a de Sitter universe or a cosmological 
constant dominated universe i.e. 

H = ^ ui = UJ2 = UJ3 = -1. (21) 

Moreover H > corresponds to an accelerating universe. This situation arises in our 
model when 

(l + cJi)Xi < 0, (22) 

(l + cu2)X2 < 0, (23) 

(1 + ^3)^3 < 0. (24) 

The above Eqs. (22), (23) and (24) yield respectively 

u,<-l, A,> ^^^^-^^^^^^~"^\ (25) 

UJ3 - UJ2 
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, / 1 \ ^ 3(^2 -^l) (^3 -^2) .^^x 

UJ2 < -1, A2 > , (ioj 

Ui - Us 

/ 1 X ^ 3(^1 -^3) (^3 -^2) 
cjs < -1, A3 > . (27) 

U2 - UJi 

Thus H will be positive when all the uJi < —1. This result corresponds to the phantom 
energy dominated universe. Note that the coupling parameters Aj have to be positive both 
for the decelerating and the accelerating universe. This result turns out to be consistent 
with [22] that coupling parameters cannot be negative to avoid violation of second law 
of thermodynamics. Moreover the same investigation shows that small positive values 
for the coupling parameters are motivated from the empirical results. Therefore the 
universe evolves from the earlier quintessence to cosmological constant and then later to 
the phantom energy dominated universe. Also the conditions (18)-(20), and the analysis 
of this section determines the super-acceleration and not the acceleration. 



2.2 Behavior of deceleration parameter 

The deceleration parameter g ~ if + H"^ = —SnG J2l=i[p'i{l + The parameter q 

is significant as its possible signature governs the dynamics of the universe. For instance 
g < represents the acceleration phase of the expanding universe. This slowing down in 
the expansion takes place when 

^2>-l \2> ^'^' .^"^^.^ (29) 

>-| A3> '^^^^ .^^^^^.^ (30) 

Thus q will be negative when all the Ui > — |. This result corresponds to the quintessence 
dominated universe. 

Note that the vanishing q in Eq. (17) will yield the following state equations 

2 

q = ui = UJ2 = UJ3 = --. (31) 

Moreover g > corresponds to an decelerating universe. This situation arises in our 
model when 

(^ + wi)Xi < 0, (32) 



3{ui 


- UJ3){UJ2 


- UJi) 




UJ3 — ^2 




S{UJ2 


- UJi){uJ3 


-UJ2) 




iOi - UJ3 






- UJ3){UJ3 


- UJ2) 


UJ2 — UJl 
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{1 + UJ2)X2 < 0, (33) 
{'^ + uj,)X, < 0. (34) 
The above Eqs. (32), (33) and (34) yield respectively 

/ 2 3(^1 -0^3) (0^2 -^l) 

3 LU^ — UJ2 

^2 3(^2 -^1) (^3 -^2) 

t^2 < A2 > , (36) 

3 ci;i — Go's 

, ^ 2 3(^1 -073) (^3 -^2) 

u;3 < — -, A3 > . [6() 

3 u;2 — ^^1 

Thus g will be negative when all the coi < — |. The strength of the interaction is deter- 
mined by the coupling parameters Aj. These coupling constants are expressed in terms 
of the equation of state parameters Ui. We here stress that the exact nature of the in- 
teraction is largely unknown i.e. the mediating particles of the interaction are not yet 
identified. Any interacting dark energy model should, in principle, be motivated from 
the particle physics or the corresponding phenomenology, however there are as yet no 
sound theoretical models which could identify the particle interactions. There are some 
arguments that the phantom-like dark energy can decay into at least one ordinary particle 
and some other phantom-like particles |23]. 



3 Stability analysis 

To perform stability analysis, we write Eqs. 

P2 + 3H{1 + UJ2)p2 
P3 + 3H{1 + LJ3)p3 



(2-5) as 

= 3ifAi(p3-P2), (38) 

= 3HX2{pi-p3), (39) 

= 3iJA3(p2-pi). (40) 



Here we have used the constraint equations on the coupling parameters by: A3 — A2 = Ai, 
Ai — A3 = A2 and A2 — Ai = A3. Since the above system is autonomous i.e. there is no 
explicit time dependent term in the equations, we can analyze the system by first finding 
its critical points and later on, checking the stability of the system about those points. A 
critical point (also called fixed point) is the one that satisfies the dynamical system (like 
Eqs. 2-5) when equated to zero. Similarly, a critical point becomes an attractor if the 
solution of the system converges to that point for large values of the time parameter. Such 
kind of analysis is generally performed to check whether the dynamical system possesses 
any solution which is stable against small perturbations [251 [26|, [271 [211 |29l [30] . 
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Further, Defining the dimensionless density parameters as 

ui = ^^n,, (41) 

Per 

U2 = — = ^^2, (42) 
Per 

Us = -^ = ^3. (43) 

Per 

Using Eqs. (41-43), we can rewrite Eqs (38-40) as 

^ = 3Mi{(1+CU2)m2+ (1+cu3)m3} + 3Ai(m3-«2), (44) 

^ = 3rt2{(l + uji)ui + (1 + u;s)us} + 3A2(iii - 1*3), (45) 

^ = 3K3{(l + ^l)^il+ (1 + c^2)m2} + 3A3(k2-hi). (46) 

Above the time derivative has been replaced with the differentiation with respect to 
a new parameter x — In a. This parameter corresponds to the number of e- foldings 
which is convenient to use for the dynamics of dark energy. By equating above three 
equations to zero, we obtain the critical points (see below). To find whether the system 
approaches to any critical point, we check the stability about these points, by producing 
small perturbations 6ui, du2 and 6U3 around the critical point {ui^,U2^,U3j i.e. 

ui = ui^ + 5ui, U2 = + 5u2, U3 = ^3^ + Sua- (47) 

Now in order to linearize our dynamical system, wc take the 6 variation of the above 
equations (44-46) and evaluate them at the critical points {{ui^.,U2^.,U3^.)}, where j is 
an indexing parameter: 



d6ui 
dx 

d5u2 
dx 

dSus 
dx 



35ui{{l + uJ2)u2^. + (1 + <^3)^i3c^ } + 2>5u2{ui^. (1 UJ2) - Ai} 
+2>5us{ui^. (1 + ujs) + Ai}, (48) 

3(5m2{(1 + ^l)uic^ + (1 + '^3)«3c^} + ^5U2,{U2^. (1 + OJ^) - X2} 

+35Mi{M2,.(l + a;i) + A2}, (49) 
35rt3{(l a;2)M2,. + (1 + } + 'i5ui{u2,^. (1 Ui) - A3} 

+35M2{M3e.(l+^2) + A3}. (50) 



We then next construct a matrix consisting of the coefficients of the perturbation pa- 
rameters and compute the eigenvalues of that matrix. The stability of the critical points 
depends on the nature of the corresponding eigenvalues and it leads to three specific cases: 
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if the real parts of all the eigenvalues are negative, then the following critical point is a 
stable node; if the real parts of the eigenvalues are all positive, then that critical point 
is an unstable node, in all other cases, the critical points will be saddle points. In our 
analysis, we shall be interested in only those critical points that yield stable nodes. 

In order to obtain some physical interpretation of our results, we assign dust (or matter), 
radiation and dark energy with subscripts 1, 2 and 3 respectively. Thus the corresponding 
EoS parameters take values cui — 0, 002 — i/^ and 003 — —1/3 or —1 depending whether 
it is quintessence or cosmological constant. From here onwards, we divide our dynamical 
system (44-46) into two parts: one for los — —1 and other for co^ — —1/3. 

It is not possible to obtain the critical points of the full system (44-46) . The only way 
out is by choosing at least one parameter Aj to be zero. In the next sections, we shall 
use only two coupling parameters while taking the third one to be zero. For instance, 
taking Ai = implies that radiation and dark energy are mutually non-interacting while 
there arc interactions between dark energy and matter and between matter and radiation. 
Similar interpretations can be made when choosing either A2 or A3 to be zero. 

3.1 Critical points for Ai = 

Note that the following critical points are obtained by setting Ai = 0. 
3.1.1 Quintessence uj^ = —1/3 

The dynamical system yields the following critical points {{ui^.,U2^.,U3^.),j = 1,2,3,4}. 



^ic, = 0, (51) 

U2^, = 1^2, (52) 
3 

^3ci = -4^3, (53) 

= ! 24A3 ( " + ^^(^^ " ^^'^^ + + 

-2 A2 ( - 35 A3 + v'l6Ai + 568A2A3 + Ai) ) , (54) 

U2^, = ^(4A2-A3 + ^16AiT568A^A^TAi), (55) 

U3^^ = ^ ( - 4A2 + A3 - Vl6Ai + 568A2 A3 + Ai) , (56) 
ui^, = 4SA -24A (-8-^2 + A3(A3 + Vl6Ai + 568A2A3 + Ai) 
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^2c3 


32 




1 


^3c3 


16 




= 0, 




= 0, 


^4c3 


= 0. 



+2A2( - 35A3 + ^IGA^ + .568A2A3 + Ai)), (57) 
^2 - A3 - ^16Ai + 568A2A3 + Ai) , (58) 



4A2 + A3 + V 16Ai + 568A2A3 + Ai) , (59) 

(60) 
(61) 
(62) 

Now we substitute above critical points in Eqs. (38-40) and obtain eigenvalues. Since 
analytical expressions for the eigenvalues are not possible, so we shall resort to estimate the 
numerical values of eigenvalues for specific choices of coupling parameters. Fixing A2 = 0.5 
and A3 = 0.7 (both positive), A2 = 0.5 and A3 = —0.7 (one positive) and A2 = —0.5 and 
A3 = —0.7 (both negative) no stable node is obtained from the above critical points. 
However, for A2 = —0.5 and A3 = 0.7, the first critical point {(^i^ ,M2c ^u^,^ )} becomes a 
stable node. 



3.1.2 Cosmological constant = — 1 



= A3, (63) 

U2., = 0, (64) 

^iSei = ^^3, (65) 

uu, = 0, (66) 

^2e, = 0, (67) 

u^., = 0. (68) 



For all the choices of A2 and A3 as adopted in the previous subsection, no critical point 
arises as a stable node. 

3.2 Critical points for A2 = 

Setting A2 = 0, we obtain the following critical points: 

3.2.1 Quintessence a;3 = — 1/3 

«ic, = -^Ai, (69) 
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U2., = 0, (70) 

«3e, - A3, (71) 

ui,^ = ^ ( - 9Ai - 2A3 + sJ^lXl + 1476Ai A3 + 4Ai) , (72) 

^2e, = 64(3Ai - 2A3) (^^-^^ - 2A3( - 2A3 + VsiAf + I476A1A3 + 4Ai) 

-3 Ai (72 A3 + ^81A? + I476A1A3 + 4Ai) ) , (73) 

= ^ (9Ai + 2A3 - sJsiXl + I476A1A3 + 4Ai) , (74) 

= ^( - 9Ai - 2A3 - ^SlAf + I476A1A3 + 4Ai), (75) 
^2.3 = 64(3Ai - 2A3) (^^-^^ + 3Ai( - 72A3 + VsiAf + I476A1A3 + 4Ai) 

+2A3 (2 A3 + ^81A? + I476A1A3 + 4Ai) ) , (76) 

«3e3 = ^(9Ai + 2 A3 + ^JslXl + I476A1A3 + 4Ai), (77) 

ui., = 0, (78) 

U2., = 0, (79) 

^3e, = 0. (80) 

For Ai = 0.5 and A3 = —0.7, the stable node arises at second critical point {{ui^^,U2^^.,U'i^^)}. 
For other values of Ai and A3, no other stable node arises. 

3.2.2 Cosmological constant a;3 = — 1 



uu, = 0, (81) 

^2e, = -^A3, (82) 

uz., = -^A3, (83) 

^ic. = 0, (84) 



U2^^ = 0, (85) 
uz^, = 0. (86) 



Here for Ai = 0.5 and A3 = —0.7, the stable node arises at second critical point {{ui , U2 , 
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3.3 Critical points for A3 = 

Now we perform similar analysis for A3 



= 0, we get the critical points 



3.3.1 Quintessence uj^, = —1/3 

= \k (87) 

U2., = -A2, (88) 

«3., = 0, (89) 

^( - 9Ai - 8A2 + ^81 Af + 4176A1A2 + 64Ai), (90) 

^ (9Ai + 8A2 - ^81Af + 4176A1A2 + 64Ai) , (91) 



ii3e, = 96Ai - I28A2 ^ 396A1A2 + 2>2\l - 3Ai^81A^ + 4176A1A2 + 64Ai 



-4A2 V81Af + 4176A1A2 + 64Aij , (92) 
: ^( - 9Ai - 8A2 - ^jsiXl + 4176A1A2 + 64Ai), (93) 
: -^(9Ai + 8A2+ ^81A2 + 4176AiA2 + 64Ai), (94) 

^3c3 = \ X (^^^1 + 396A1A2 + 32A^ + ZXJsiXl + 4176A1A2 + 64Ai 

^ 9oAi — I28A2 ^ 



+A\2\jsi\l + 4176A1A2 + 64Ai), (95) 
uu, = 0, (96) 
«2., = 0, (97) 

Ii3e, = 0. (98) 

Among the above critical points, the only stable node is produced for the first one for 
values Ai = —0.5 and A2 — 0.7. 



3.3.2 Cosmological constant 0^3 = —1 



4^- 
-A2, 

0, 



(99) 

(100) 
(101) 
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7A iA2(-9Ai + I 6A2; 
4(3Ai +4A2) 

«i.3 = 0, (105) 
= 0, (106) 
= 0. (107) 



Among the above critical points, the only stable node is produced for the first one for 
values Ai = —0.5 and A2 = 0.7. 

In figures 1-6, we provide the schematic representation of triple interacting fluids, taking 
two at a time, with various choices of the coupling parameters and initial conditions. The 
initial conditions are chosen arbitrarily but these must satisfy the constraint mi(0) + 
^2(0) + ^3(0) = 1 (the Friedmann equation). In figures 1 and 4 , we assume first and 
second interacting fluids whereas in figures 2 and 5 and 3 and 6, we assume second and 
third and first and third interacting fluids respectively. The three solutions ui, U2 and 
always satisfy the constraint mi(0) + ^2(0) +^3(0) = 1. 



4 Conclusion and discussion 



In this paper, we have attempted to resolve the cosmic triple coincidence problem which 
naively asks why the energy densities of the three major ingredients of the cosmic com- 
position namely matter, radiation and dark energy are of same order at current time. 
Rephrasing, why uj has evolved to —1 in recent times. We here point out that the EoS 
p = ujp used in this paper is not, in general, a true EoS for any cosmic fluid. Rather it 
is a phenomenological relationship suitable for the configuration. The actual EoS may 
not be that simple and may have dependencies on various other cosmological parameters 
like redshift, time, Hubble parameter and its derivatives etc [8]. However to a first order 
approximation, it may be taken for such analysis. Further, the standard non-interacting 
FRW model cannot resolve the coincidence problem since it predicts a hierarchial system 
in which radiation density decreases faster compared to matter, while density of dark en- 
ergy remains either constant (if it is cosmological constant) or increases (if it is phantom 
energy). Therefore it contradicts with the observations where all the three components 
have equivalent densities. Model of interacting dark energy, which is a modification of 
the FRW model, has enormous potential to explain this cosmological conundrum. It nat- 
urally predicts that if the cosmic fluids interact with each other, it leads to a scenario 
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compatible with the observations [23] . 

We have found that in the interacting fluid model, the three fluids can achieve super- 
negative equation of state. If matter and radiation are among the components then they 
will induce negative pressure along with the dark energy to produce accelerated expansion. 
This result has earlier been shown for the interacting Chaplygin gas model in [7j. 
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Figure 1: The phase diagram of the triple interacting fluid model with the choice of couphngs 
[Ai = 0.2, A2 = 0.04, A3 = 0.06], X = — 1,..,1 and EoS parameters [ loi = 0,uj2 = 1/3, uj^ = 
— 1]. Chosen scene is between ui and U2. The curves correspond to the initial conditions 
ui{0) = 0.2,?X2(0) = 0.6,^3(0) = 0.2 (brown); ni(0) = 0.4,^2(0) = 0.2,M3(0) = 0.4 (red); 
ui{0) = 0.3,W2(0) = 0.4,7/3(0) = 0.3 (green). 
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Figure 2: The phase diagram of the triple interacting fluid model with the choice of couplings 
[Ai = 0.2, A2 = 0.04, A3 = 0.06], X = — and EoS parameters [ wi = 0,u;2 = 1/3, = — 1]. 
Chosen scene is between U2 and U3. The curves correspond to the initial conditions as taken in 
Fig.l. 
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Figure 3: The phase diagram of the triple interacting fluid model with the choice of couphngs 
[Ai = 0.2, A2 = 0.4, A3 = 0.6], X = — 1,..,1 and EoS parameters [ wi = 0,u;2 = 1/3, = — 1]. 
Chosen scene is between ui and U3. The curves correspond to the initial conditions as taken in 
Fig.l. 
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Figure 4: The phase diagram of the triple interacting fluid model with the choice of couphngs 
[Ai = 0, A2 = 0.4, A3 = —0.4], X = — 1,..,1 and EoS parameters [ toi = 0,uj2 = 1/3, = —1/3]. 
Chosen scene is between ui and U2- The curves correspond to the initial conditions as taken in 
Fig.l. 
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Figure 5: The phase diagram of the triple interacting fluid model with the choice of couplings 
[Ai = 0, A2 = 0.4, A3 = —0.4], X = — and EoS parameters [ toi = 0,uj2 = 1/3, = —1/3]. 
Chosen scene is between M3 and U2- The curves correspond to the initial conditions as taken in 
Fig.l. 
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Figure 6: The phase diagram of the triple interacting fluid model with the choice of couplings 
[Ai = 0, A2 = 0.4, A3 = —0.4], X = — 1,..,1 and EoS parameters [ toi = 0,uj2 = 1/3, = —1/3]. 
Chosen scene is between ui and U3. The curves correspond to the initial conditions as taken in 
Fig.l. 
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